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Abstract 

It is demonstrated that the planar reorientation of 
a free-free beam in zero gravity space can he ac- 
complished by periodically changing the shape of the 
beam using internal actuators. A control scheme is 
proposed in which electromechanical actuators excite 
the flexible motion of the beam so that it rotates in 
the desired manner with respect to a fixed inertial 
reference The results can be viewed as an extension 
of previous work to a distributed parameter case. 

1. Introduction 

Following [8], we introduce the concept of a de- 
formable body, for which distances between the points 
of the body can change during the motion. Lx- 
amples of deformable bodies include both ''imped 
and distributed parameter systems such as multilmk 
rigid body interconnections and structures with dis- 
tributed flexibility. The orientation of a deformable 
body with respect to a fixed inertial reference can be 
specified by a choice of body frame In general, thexe 
are many ways to choose a body frame. For exam 
pie. in the case of planar motion a body frame ^can be 
identified with any two distinct points in the 
The shape of a deformable body can be specified in 
terms of the position of the body relative to the body 
frame Thus" an arbitrary motion of a deformable 
body can be separated into rigid body motion and 
shape change. 

Assume that both linear and angular momenta about 
the center of mass of the body are conserved and 
eoual to zero. These conditions hold if the body is 
In a circular orbit around the Earth or is in a free 
fall As a consequence of angular momentum conser 
vation shape change and the rigid body motion are 
coupled. This coupling is inherently nonlinear In 
., . .L.u, iptere sted in induc ing a ro^ 
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In this paper we extend the aforementioned reorien- 
tation strategies to the case of flexible bodies. In 
particular we are interested in a planar reorientation 
of a free-free beam in space using only electrome- 
chanical actuators. These electromechanical actua- 
tors e g piezoelectric or shape memory actuators, 
do not change the angular momentum of the tree 
free beam but. can be used to change the shape ot the 
beam in a periodic way. Assuming that the .angular 
momentum of the beam is always zero oscillations 
in the shape of the beam can cause a rotation ot the 
beam with respect to a fixed lnert *f l , refe ^ c ® nl 
rotation of the beam over one period depends °nly on 
the shape of the beam over one period and does not 
depend on the length of the period; hence ^{f 1 Pjf 
nomenon is referred to as a geometric phase change. 

The extension of existing strategies to the free-free 
beam case is not straightforward for sevejalreasons^ 
Classical models of uniform free-free flexible beamb 
in zero gravity space result in complete decoupling 
of rigid body motion and flexible motion. Higher 

ordef nonlinear coupling between rigid inexact 
and flexible motion is captured in geometrically exact 
beam theories [9]. The resulting models, however, are 
complicated. The free-free beam is an infinite dimen- 
sional superarticulated system. Thus, an arbitrar 
shape change cannot be produced by a finite number 
of aerators 8 In addition "the body frame of the beam 
needs to be chosen so that the shape change is inde- 
pendent of the rigid body motion. Such a choice of 
Pody frame is natural for lumped parameter systems 
since variables specifying orientation are ignorable 

In this paper, we first address basic modeling is- 
sues" ^"dynamics which determine the shape o 
the free-free beam are assumed to be characterized 
bv the Euler-Bernoulli equation, including material 
damping, with appropriate boundary conditions. The 
higher order coupling between the rigid body motion 
arfd the flexible motion is captured using the angular 
momentum expression which includes rotatory iner- 
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frame with its origin fixed at *he origin of the iner- 
tial frame such that the vectors (i. k) lie in the plane 
( e i , e 3 ) and j = e? • The straight line passing through 
the origin in the direction of vector k is called the 
reference line. Let the beam initially be at rest in a 
straight line configuration aligned with the reference 
line. Then, the location of each point on the line of 
mass centroids of the beam can be described in terms 
of the parameter s € [— L, L\. This parameter .s can 
be viewed as a label for each of the orossections. We 
assume that as the beam deforms the shape and the 
area of the crossections remain invariant. Following 
other researchers [1. 6, 9] we introduce three func- 
tions • [— £. L] x ?Ji — ft and rf's./) : 

[—£,£] x 'ft — T x such that (u(.s. /) + s. y(s. t)) de- 
fine the coordinates of the line of centroids in the 
deformed configuration with respect to the moving 
frame (i,£) at time t. The angle t*(s.n between the 
normal to the crossection at s and specifies the 
orientation of the crossection. The normal to the 
crossection at s is denoted by £3. We define the mate- 
rial basis (t 1 T2T3) to be orthonormal so that t\ lies in 
the plane ( e 1 , e 3 ) - The crossection itself can be asso- 
ciated with the set of points (£1,^2) in a compact set 
A C S 2 such that +£2*2 +{u(s, t) + s)k + (y(s, t))i 

gives the location of any point on the beam as £1 and 
£2 vary through .4 and s varies from — L to L. 

Since the origin of the inertial frame is fixed at the 
center of mass of the beam we obtain 


separation between the motion which determines the 
shape ot the beam, given by y(s.t). -L <s < L. and 
the rotation of the beam as a whole, given by <Ht). 



Fig 1. Planar Beam Model 


We next develop a kinematically exact expression for 
the angular momentum of the free-free beam. Let 
£1 1 #T) be the vector from the origin of the 

inertial frame to a point ($,£1^2) on the beam at 
time t; then 



p = (s sin 9 4* y cos 9 + £1 cos 0 + u sin 9 )e\ + 

(£2)^2 + (5 cos 9 — £1 sin xp — y sin 9 + u cos 0)e 3 (7) 



where 9 — 9 {t),y = y(sj) and rp = rp(s,t)- The 
(2) angular momentum about the origin of the inertial 

frame at time t is zero so that 


Let p denote the constant mass density per unit vol- 
ume of the beam. We assume that the beam has a 
symmetric crossection so that the first moment of in- 
ertia of the crossection about the line of centroids is 

J r Ki^irf6=0. (3) 



dp 

P¥ x = 0. 


{*) 


Substituting equation (7) into equation (8) and using 
equations (4) and (5) we can express 9 in terms of 
y . a and a as 


The second moment of inertia of the crossection about 
the line of centroids, referred to as the rotatory iner- 
tia, is 

(4) 

and assumed to be positive. The mass per unit length 
of the crossection is given by 


m 0 



(5) 


We define the angle 9 (t) between e 3 and k so that 
y(s,t) measured from the reference line satisfies the 
following orthogonality condition 



sy(s, t)ds = 0. 


( 6 ) 


The existence of the angle 9 (t) follows from the ge- 
ometry indicated in Fig. 1. This definition provides a 


Q _ + ha + n»o(f*« - 

jf L {-m 0 s 2 - m 0 y 2 - I 2 }ds 

where a — </> — 9 is the angle between the normal / 5 
to the crossection at s and the reference line. 

Assume that the beam is unshearable and inextensi- 
ble and that the deformations are small. This implies 
using equation (2), that 

u(s y t) — 0. f I'M 


and that 


assy,. 
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We use the Euler-Bernoulli beam model to charac- 
terize the shape of the beam [3]. Thus y{s t t) satisfies 
the Euler-Bernoulli equation of the form 


m 

moytt + iytssss + EIyssss = 6 '(s-Sj) ill’ 

j = i 


I 



with the boundary conditions • 


y,A~L) = Vss(L) - O' 

t/j j * ( L) = J/aJj(^) = 0 

where I = h/p< E is Young’s elasticity modulus. 6' is 
the distributional derivative of the delta function and 
where for simplicity we assume Kelvin- Voigt damp- 
ing with a positive damping coefficient 7 . In addi- 
tion, y(s,t) must satisfy conditions (1) and ( 6 ). In- 
ternal bending torques fj(<), j = 1 m are pro- 

duced by m point actuators located at s = Sj on the 
beam where sj € [-L. L). These actuators change the 
shape of the beam but at the same time preserve the 
angular momentum. Although such actuators are ca- 
pable of inducing relatively small displacements one 
can excite the beam periodically at a frequency near 
one of the lower resonant frequencies of the beam to 
obtain relatively large periodic shape change. 


(13) 

( 14 ) 


Using expressions ( 6 ), (10) and ( 11 ) in equation (9) 
we obtain . 


r + f^ L may 2 ds 

where t = §m 0 L 3 + 2/ 2 L. This expression demon- 
strates the nonlinear coupling between the beam s 
shape and its rigid body motion. Expression (15) is 
non-integrable in the sense that if y(s,<) is a periodic 
function of time, the integral of 9 over one period is, 
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3. Asymptotic Reorientation Maneuvers 

The goal is to accomplish asymptotic maneuvers, i.e. 
starting with 0(*o) = #0i y( s t o) = 0 ) = 0 we 

want to rotate the beam so that 9(t) — $ 4m y(sj) — 0 
and t/t(M) — 0 as f — * x. ^ or some desired angle 9 {i . 

Consider the periodic excitation of the beam at a sin- 
gle frequency as 

Vj{t) - ij + Vj cos (ut),j = 1,2 rn ( 19) 

Since the shape space dynamics of the free-free beam 
is asymptotically stable, the steady-state motion of 
the beam is given by 

q,{t) = lx + n, cos(w< + ©i ) ( 29) 

where the parameters t t , a; and d>, ran be expressed in 
terms of vj and t ; ° according to equation (17). The 
excitation function (19) should be sufficiently small 
so that the Euler-Bernoulli model for the shape space 
dynamics remains valid. Substituting equation (20) 
into equation ( 18) and integrating over one period we 
obtain the steady-state change in angle 9 over one 
period is given by 

f <o cos(u^f + Xo )dt (21) 

J 0 1 + m(cosM + Xl) + <2 cos(2w< + Xi) 

for constants coi <u> XOi Xl and Y 2 Expres 
sion (21) implies that, in g ene ral, the change lnjmgle 
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non-straight line reference configuration. It follows 
from expression ( 18) that in order to rotate the beam 
in the opposite direction it is sufficient to reverse the 
signs of v* and Vj . 

We are now in a position to formulate a specific con- 
trol strategy to accomplish the desired asymptotic 
maneuver. Starting at rest with 0(fo) - 0o applica- 
tion of control law (19) results in a nonzero geometric 
phase change over one period. By repetition of cycles 
of motion as many times as necessary the beam can 
be caused to rotate closer and closer to 0 As 0(t) 
approach 0 j we can reduce the amplitude of the oscil- 
lations to zero in a way so that 0 (/) — 0 * as t — x 


The proposed control law is of the form 

t j{t) = Sic [t?° 4- ij cosM)] J = 1 n} (23) 

where 2 { k ~ l ) w P < t — f 0 < f. — 1 . *J. . . that 

u/ — u j 

is, the control excitation is an amplitude modulated 
function, where t'°, v* are constants and s k denotes 
the scalar amplitude modulation sequence that de- 
fines the control excitation on the A:-th cycle. Each 
cycle is exactly p periods. The constants w. Oj , vj can 
be chosen nearly arbitrary, although one approach is 
to choose Vj, v* to maximize geometric phase expres- 
sion (22) where a, „ / t , <pi, i = 1 , . . . , are related to v'j , 
Vj 1 . j — l,...,m according to expressions ( 20 ) and 
(17), and v*j , v* are constrained in norm. In terms 
of Vj, vj . j =. 1 this is a constrained mathe- 

matical programming problem which is linear in Vj 
(for fixed t>") and quadratic in v* (for fixed vj). We 
will subsequently denote the maximum value of this 
constrained optimization problem as A0T 


The modulation sequence s k + \ is defined in terms of 
an average of 0(t ). over the Ar-th cycle, that is 




1 
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(max 0(f) 4- min 0(0) 


(24) 


where the maximum and minimum are over 

-< fc-li ^ -H*t p 


< t - tn < 


We also introduce two 


auxilary variables — 0 r) and Sq = sign ( ) . 
We express s k in terms of 0 k v J { an d -fc-i ^ indicated 
below: 


(Al) Compute 


n e = 



(A2) In case \r k \ > if r k and s^-i have 

the same signs then Sk = |f t- 1 jsign(rjt); if 
r k and s k -i have opposite signs then s k = 
7 i |cib_ i |sign(r* ), where 0 < 7 i < l- 

(A3) If 0 < \r k \ < \s k -\\ then s k = 72 r*, where 

0 < 72 < 1. 

(A4) If r k = 0 then s k = 


Proposition 4.2 If the proposed control law ls 
of the form ( 2 :j) where -: k is selected according to 
steps (A1)-(A4). then 

lim 0 k ve = 0,i. lirn e k — (j. 

v -x k — x. 


Sketch of the Proof. The sequence \-; k j is nun - 
increasing and bounded on [ 0 , l]. Therefore, there 
exists b 6 [ 0 , 1 ] such that b = inf ^ I - It - an !>«- 
shown that by construction of the sequence b mii>t 
be zero. 


Since \~ k \ — 0 then q t (t) — 0 and q t — 0 as t — x 
By continuity 0(0 — 0 con for some constant if' n a 
t — x It can be shown that 0 con = 0 ti . _ 


Finally, it follows from equations (24) and ( 20 ) that 

lim 0(0 = 0*, lim ( ) = 0, -L < s < L 

The controller which we have constructed has two 
functions. Its main function is to excite the oscil- 
lations of the beam in such a way that the beam 
rotates in the desired sense. Subsequently, the con- 
troller serves to suppress the vibrations previously ex- 
cited so that the free-free beam comes to rest with a 
desired orientation. Note that control law (23) is a 
non-smooth feedback control law [ 2 ]. 


4. Numerical Example 

Consider a beam with half-length L = l[m]. den- 
sity per unit volume p = 1400[A:^/m 3 ] and square 
crossection with the side size R = 0.1 [ 771 ]. Young s 
modulus of the beam is E = 3.0 x 10 6 [iV/m 2 ] and 
the Kelvin- Voigt damping coefficient is 7 = 0.2. Tw>» 
actuators are installed near both ends of the beam 
at = -0.9[m] and s 2 = 0.9[m]. The maximal 
torque each of the actuators can produce is equal . 
lOU[.Ym]. The excitation frequency w = 13[/7_- is 
lected to lie between the first 10 .6[H z] and the second 
29[H z] resonant frequencies of the beam; v® and >:'! 
j — 1.2 are chosen using expression (22) to maximize 
the geometric phase change over one period. For thi> 
example we choose p = 0 and y\ = 72 = 0.9. Tie- 
first four elastic modes of the beam are used in «»ur 
simulation. 


We want to rotate the beam from 0q = O.lfrrn/] at 
t = 0[sec] to 0<i =: 0[ra<f]. The dependence of tli^ 
angle 0(£)[rarf] on time t[sec] is shown for a part - f 
the maneuver in Fig. 2. In this case the georm*tri< 
phase change over one period in steady-state pre- 
dicted by expression (22) is equal to -2.7465 x lU“ + 
[rad] whereas its actual simulation value is equal r 
—3.0411 x 10“ 4 [rad\. The dependence of the m<, ;- 
ulation parameter c on time is shown in Figure ^ 
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Fig 2. Asymptotic Reorientation Maneuver 


beam by internal actuators are derived. Finally, a 
control strategy for a planar asymptotic reorientation 
maneuver is developed. 


A general treatment of the interplay between defor- 
mations and rotations of deformable bodies is given 
by Shapere and Wilczek [8]. Reyhanoglu and Mc- 
Clamroch [7] have developed a framework for reorien- 
tation of mult.ibody systems in space. In this paper, 
we have used the framework developed by Shapere 
and Wilczek for the specific problem of reorientation 
of a free— free beam in space; our results represent, in 
a certain sense, the limiting case of the multibody re- 
sults obtained by Reyhanoglu and McC lamroch when 
tha nnmhpr of bodies increases without limit. 


\|though our study in this paper has been concerned 
with the ideal case of reorientation of a free-free beam 
in space we note that the same ideas are applicable 
to reorientation of a wide class of deformable space 
structures, using only actuators embedded into the 
structure. In this sense, smart structures technology 
can be used to accomplish a variety of efficient reori- 
entation maneuvers for space structures. 
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Fig 3. Amplitude Modulation Sequence 


5. Conclusion 

In this paper the angular momentum expression for a 
planar free-free beam in space is derived. It is shown 
how the general motion of the beam can be separated 
into rigid and elastic motions. The change of shape 
ptibe beam is described by the E uler-Bemoull i equa- 
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